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Background from the WMAP Temperature Maps 
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Abstract 

In this paper, we study how to predict the polarization of the Cosmic Microwave 
Background using knowledge of only the temperature (intensity) and the cross-correlation 
between temperature and polarization. We derive a "Wiener prediction" method and 
apply it to the WMAP all-sky CMB temperature maps and to the MAXIMA field. 

1 Introduction 

Fluctuations in the temperature of the Cosmic Microwave Background (CMB) have been 
used to map the surface of last scattering for the last decade , since the first results from the 
DMR instrum ent on the COBE Sat ellite rtSmoot et al. Now, we have the first results 

from WMAP (|Bennett et al.l . 1200,1 references therein and elsewhere in this volume) which 
confirm the CMB anisotropy results gathered over the intervening decade. But the WMAP 
data promise more than just co nfirmation: they als o offer the first high- sensitivity analysis 



of the polarization of the CMB (jKogut et al.L 120031 ) [although PAST ca n justly be credited 



with the first detection of cosmological polarization (jKovac et al. ]. l2002h ]. 

The polarized CMB provides independent information on cosmological physics; it traces 
the flow of the plasma at the surface of last scattering. This flow itself can be decom- 
posed into an irrotational component, corresponding to the action of gravity by the den- 
sity perturbations responsible for the temperature fluctuations, and a rotational compo- 
nent which can be produced by a background of gravitational waves {e.g., from inflation). 
These comp onents result in patterns o f CMB polarization known, respectively, as "grad" 
and "curl" (iKamionkowski et al. . 1997 ) or, in analogy with electrodynamics, "E" and "B" 
( Zaldarriaga and Seliakl . Il997l ). The latter (curl or B) thus has the promise of exploring the 
physics of the early Universe, but the amplitude of the signal remains well below current 
detections. Here, then, we concern ourselves with grad or E modes. 
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However, the WMAP team has yet to produce actual maps of the polarized component of 
the CMB. Here we investigate a statistical technique to predict the polarization component 
based on knowledge of the auto- and cross-correlations of the components. 



2 Formalism 

In this section, we ask the question: if we know the temperature pattern on the sky, as 
well as the cross-correlation between the temperature and the polarization, what is our best 
guess at the latter? There are many ways to tackle this question, mathematically if not 
philosophically equivalent. Here, we take a Bayesian approach, and so start with Bayes' 
theorem, 

where 9 labels the parameters we are trying to determine, D the data, and I the back- 
ground information we are bringing to the problem. P{a\b) is the probability (or probability 
density) of a given b, so P{6\I) is the "prior" for 9 and P{D\9I) is the likelihood — the 
probability of getting the specific data D given a fixed set of parameters. The left hand side 
is the posterior probability, and the denominator just enforces the normalization condition: 
P{D\I) = Jd9P{9\I)P{D\9I). 

How do we apply this to the CMB problem? Let us start with a model for the data, Ap 
as a sum of signal, Sp and noise. Hp, 

Ap = Sp + Up , (2) 

where Ap gives the value of the data labeled by p, which can include pixel number, frequency 
of observation, polarization state, etc. We will take the noise to be described by a zero-mean 
Gaussian distribution with 

(up) = 0; {upUp') = Nppi (3) 

where angle brackets are averages over the distributi on of the data, and N^w i s the noise 
covariance matrix (which we assume is known, but see lFerreira and Jaffd ( 20001 ) : iDore et al. 



( 200l[ )). We will also take the signal to be a zero- mean Gaussian, with analogous properties 



to those above, and a signal covariance given by {spSpi) = Spp'. 
Thus the actual likelihood function is given by 



P(A|,SiV) = ^exp 



-^A\S + N)-'A 



(4) 



where we use matrix notation, f refers to the Hermitian conjugate, and vertical bars to the 
determinant. We can then use Bayes' theorem to ask for the posterior distribution of any 
quantity related to the signal, Sp. 

The most obvious quantity we might estimate is the signal itself, 

P{s\ASN) oc P{s\SN)P{A\SN) 
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oc exp[-^s^S-h-^A^{S + N)-^A] 



oc exp 



{s\A)ya 



t r^-i 



(^|A)) 



(5) 



where the second hne comes from assigning a Gaussian prior with variance Spp' to the signal, 
and the final line comes from completing the square, and we have defined 



{s\A) = S{S + N)-^A and Cs = {Ss 6s^A) = S - S{S + N)'^ S 



(6) 



which give the Wiener filter itself, and the variance about the mean 



The W iener filter 



Bunn et al 



but 



has been used in the context of the CMB many times before {e.g. 
primarily as a method of smoothing or for the separation of various foreground components 
in the signal. 

But we can ask a somewhat more general question: what is the distribution of some 
quantity, for which the prior information is that {epCpi) = Eppi and that it is correlated 
with the signal, {cpSpi) = Xppil Again, the result is a Gaussian distribution, this time with 
mean and variance given by 



{e\A) = X{S + N)-^A and = E - X{S + N)-^ X 



(7) 



In fact, all of these results can be compressed into two equations, using the notation of 
baitoL^ B. already used to some extent above. 



{x\d) 
{Sx 6x^\d) 



{xd''){dd^) , and 
(xx^) - {xd^){dd^)-\dx^) 



where d = s + n and x is correlated with s, so (xd^) = (xs^), which is the cross-correlation, 
X, above. If x = s, this reduces the original Wiener filter, Eq. IHl above. 

Information about the variance of the Wiener filter is crucial. The Wiener mean itself, 
{x\d), is over-smoothed — low in regimes of low signal-to-noise ratio. A more "typical" signal 
value is provided by an actual realization from the full posterior distribution: the mean plus 
a realization of the correlation matrix. 



3 Applications 
3.1 Full-sky maps 

The action of the Wiener filter prediction is clearest when we consider its action on an all-sky 
map with uniform (or axisymmetric) noise. We start with a noisy temperature map, which 
we transform into its spherical harmonic components, dim = aim + nim, where aim gives the 
CMB signal and n£m the noise. The Wiener filter and variance are given by 
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{Eem\dem) — — Ardgm 



r<TE 

2 



{SaemSaem\dem) — Ce — 2^ 



Bid + m 

{^Eirn^Eirn\dlm) = Cg 



B]Q + N, 
riEEn 



{Saem5Eim\dem) — — — — — - , (9) 

where we have introduced the experimental beam, Bi, the temperature power spectrum, 
Ci = (la^mP), the El-mode polarization power spectrum, Cf = (jE'^mP) and the E-mode 
polarization/temperature cross power spectrum, Cj^ = {al^Eim)- These power spectra 
are diagonal in i, and we impose a similar constraint on the noise power, which we denote 
Ne = (|?i^mP)- Note that this is an unphysical condition, corresponding to isotropic noise 
over the sky; if we restrict ourselves to high signal-to-noise experiments and large scales. In 
these equations, we show the Wiener filter for the temperature itself and for the polarization, 
and, by extension of the equations from the previous section, the covariance between the 
two. 

We note in passing that the formalism also correctly accounts for the statistics of the 
polarization "B-modes". Since the cross-spectrum Cj^ = 0, we just have {af^\dim) = with 
{dBirn^Birnldim) = Cf — the data give no new information about them. 

We show these formulae for the Wiener filter of the multipole coefficients themselves, but 
in fact we will always want to smooth the resulting map with some beam (although it is not 
necessary to use the same as the original experimental beam, Bi). 

We can apply these formulae directly to the temperature data from the WMAP satel- 
lite. Our requirement of isotropic noise should be adequate for I < 400 where the noise is 
ne gligible. In fact , in th e following we use the Wiener-filtered temperature maps provided 
by iTegmark et al. Note that starting from the temperature Wiener filter enables a 



slight simplification: 

(jTE 

{Elm\dlm) = {0'im\dlm) ■ (10) 

Starting from these Wiener-filtered data, which are in fact available in pixel space, we 
apply these formulae by going back and forth between pixels and multipole coefficients 
by performing a spherical-harmonic transform using routines from the HEALPix package 



fjGorski et al.l . Il999).^ In Figure ^ we show the Wiener-filtered map of the WMAP data. 



with the polarization superposed on the temperature. 

In Figure El we show 4 realizations of the temperature and polarization drawn from 
the variances of Eq. IHl It is evident that, in fact, the polarization field contains a large 
uncorrelated component. We can understand this by examining the correlation coefficient. 



^http://www.eso.org/science/healpix/ 



- M!)! (11) 
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Vpi'oj/™AP/tits/WTE_map_0THTfiener_256_64.fits: Temperature + Polarisation 
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Figure 1: Wiener-fil tered map of th e CMB temperature and polarization using the WMAP 
data as prepared by Tegmark et al. ( 20031 ). 
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Figure 2: Two realizations of the map of the CMB temperature and polarization using the 
Wiener-filtered map Figure Q 



5 



This correlation coefficient gives the fraction of the power at a given multipole that is corre- 
lated between the 13-mode polarization and the temperature. In the signal-dominated regime 
(here, low tj this gives the correlated fraction of the power in the Wiener filter, as well. In 
addition, the quantity 1 — rj gives the fractional amplitude of the variance of the Wiener 
filter. Thus, if ~ 1, a realization of the Wiener variance is small compared to the Wiener 
mean, and the mean is a good indication of the actual signal, but if is closer to 0, dif- 
ferent realizations can look very dissimilar. We plot the correlation coefficient in Figure El 
We see that the correlated fraction depends strongly on i (and indeed drops to zero where 
Cj^ = 0) and is everywhere less than about 0.8. In practice, this means that the realizations 
do vary considerably, but this can in principle be ameliorated by careful filtering of the map 
to remove those modes where CT^ ~ 0. 



3.2 Small areas 

We can also use the same technique to predict the polarization for small areas of sky, where 
the available data may be more complicated than that from WMAP. In Figure |3] we show 
a prediction of the polarization for the field probed by the MAXIMA experiment {e.g., 



Lee et al.L l200l[ ). Here, we take the full noise and signal covariance into account, and per- 
form the calculation in pixel space. Instead of performing the calculation i-hj-i, we must 
now perform full matrix operations. For the isotropic CMB temperature field with power 
spectrum Cg, 

Of A- I 

Spp' = E ^^C7,S2p,(cos V) (12) 

where Bi gives the spherical harmonic transform of the experimental beam and 9pp/ is the 
angle between pixels p and p'. Analogous formulae hold for the covariances of the polariza- 
tion components, and slightly more complicated ones for the cross-correlations amongst and 
between different polarization and temperature components. 



4 Discussion 



As mentioned above, the correlation coefficient, ri gives an indication of the degree to which 
the the polarization signal is completely predicted by the temperature signal (and the power 
spectra Cf^ ). The fact that |r^| < 0.8 is an indication that a significant fraction of the 
E-mode polarization signal is not correlated with the temperature. This leads, finally, to 
a question about the physics of the temperature-polarization correlations: the correlated 
component of the polarization and temperature is due in large part t o correlations between 
density and velocity perturbations on the surface of last scattering dK amionkows ki et al. 



19971 : IZaldarriaga and Seliakl . Il997l : I.Taffe et all . hOOd IHu and Dodelsonn2002l . Having only 



temperature data as discussed here means that the results assume such a relationship be- 
tween the two components. However, with the availability of actual polarization maps, we 
should be able to use only minimal theoretical input and yet still recover a two-dimensional 
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Figure 3: Top Panel: the correlation coefficient, r, between temperature and E-mode CMB 
Fluctuations (bottom curve); the top curve shows 1 — r^. The bottom three panels show the 
individual power spectra, Cf^ . The vertical lines give the locations where Cj^ = 0. Note 
that the model we use here has optical depth, r = 0.17, responsible for the upturn in the 
spectra and correlation at low i. 
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Figure 4: Wiener predicted polarization for the MAXIMA-1 field. 
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snapshot of the plasma flow on the surface of last scattering. We leave the details of this to 
a future work and, we hope, future data. 
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